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Abstract 

Assessing the stability of higher-dimensional rotating black holes requires a study of linearized 
gravitational perturbations around such backgrounds. We study perturbations of Myers-Perry 
black holes with equal angular momenta in an odd number of dimensions (greater than five), 
allowing for a cosmological constant. We find a class of perturbations for which the equations 
of motion reduce to a single radial equation. In the asymptotically flat case we find no evidence 
of any instability. In the asymptotically anti-de Sitter case, we demonstrate the existence of a 
superradiant instability that sets in precisely when the angular velocity of the black hole exceeds 
the speed of light from the point of view of the conformal boundary. We suggest that the endpoint 
of the instability may be a stationary, nonaxisymmetric black hole. 



1 Introduction 

Exact solutions describing higher-dimensional rotating black holes have been known for a long time 
^ but the question of their classical stability is still unresolved. There are arguments suggesting 
that a Myers-Perry (MP) black hole will be unstable for sufficiently large angular momentum in 
five |2] and higher ^ dimensions. However, a convincing demonstration of this requires a study of 
linearized gravitational fluctuations around such backgrounds. This has only been done in the limit 
of vanishing angular momentum, i.e., for higher-dimensional Schwarzschild black holes [lJ|Slin]. 

Another context in which stability of higher dimensional rotating black holes has been discussed is 
the AdS/CFT correspondence [7j. The MP solutions have been generalized to include a cosmological 
constant |H1 El • There is a qualitative argument that rotating, asymptotically AdS black holes might 
exhibit a superradiant instability ^]. The idea (inspired by the corresponding instability of a Kerr 
black hole in the presence of a massive scalar field |1 HI121[T3] ) is that superradiant perturbations are 
trapped by the AdS potential barrier at infinity and reflected towards the black hole where they get 
amplifled and the process repeats. It can be proved that no such instability is present for black holes 
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rotating at, or slower than, the speed of hght relative to the conformal boundary ^01 (essentially 
because superradiant modes do not fit into the AdS "box" IT1|) but an instability may well be 
present for more rapidly rotating holes. The only way to find out is to study perturbations of such 
black holes. This has been done for scalar field perturbations of small four-dimensional Kerr-AdS 
black holes p!5| but not for gravitational perturbations, large black holes, or higher dimensions. 

The main goal of this paper is to study gravitational perturbations of higher-dimensional rotating 
black holes. In four dimensions, gravitational perturbations of rotating black holes can be studied 
analytically. The tractability of the problem arises from two miraculous properties of the Kerr metric. 
Firstly, the equations of motion for gravitational perturbations can be decoupled and reduced to a 
PDE for a single scalar quantity J^. Secondly, this equation can be reduced to ODEs governing the 
radial and angular behaviour by separation of variables. It is known that this separability property 
of the Kerr metric does extend to some of the higher-dimensional MP metrics jl7| I19j , and MP 
metrics with a cosmological constant POj. This makes the study of scalar field perturbations in such 
backgrounds tractable. However, so far no-one has succeeded in decoupling the equations of motion 
governing gravitational perturbations of MP black holes. 

In this paper, we will make progress with this problem by considering the subclass of MP black 
holes for which the number of space-time dimensions is odd and the angular momenta are all equal, 
allowing for a cosmological constant. Such black holes are cohomogeneity-1: the metric depends only 
on a radial coordinate. In D = 2A^-|-3 dimensions, the isometry group is enhanced from R x U{1)^~^^ 
to Rx U{N+1) where R denotes time translations. The horizon is a homogeneously squashed 5^^"*"^ 
viewed as a bundle over CP^ . The symmetry enhancement allows us to classify gravitational 
pertubations into scalar, vector and tensor types according to how they transform under isometrics of 
(jpN Tensor perturbations do not exist for = 1 so we restrict ourselves to A^ > 1, i.e., spacetime 
dimensionality seven or greater. In general, the different types of perturbation are coupled together 
but we shall show that there is a subset of tensor perturbations which decouples. After separation 
of variables, the equations of motion for such perturbations reduce to a single ODE governing the 
radial behaviour. For completeness, we also consider massive scalar field perturbations of these black 
holes (for A^ > 1). We shall present a unified form for the radial equation that applies both to scalar 
fields and to the gravitational perturbations just discussed. 

Asymptotically flat, cohomogeneity-1 Myers-Perry black holes exhibit an upper bound on their 
angular momentum for a given mass. Solutions saturating this bound have a regular but degenerate 
horizon. This means that the black holes considered here behave rather differently from ones for 
which some of the angular momenta vanish, which are the ones expected to exhibit a gravitational 
instability [SUB]- So there is no a priori reason to expect any instability to exist for the MP black 
holes considered in this paper and indeed we find no instability. 

In the asymptotically anti-de Sitter case, there is also an upper bound on angular momentum 
for given mass and black holes saturating this bound have a regular but degenerate horizon. These 
extremal solutions always rotate faster than light and can be arbitrarily large compared with the 
AdS radius. The argument of 10^ suggests that black holes near to extremality might be unstable to 
losing energy and angular momentum into superradiant gravitational and scalar field perturbations. 
We shall demonstrate that this is indeed the case. Moreover, we shall show that this instability 
appears as soon as the angular velocity of the hole exceeds the speed of light, i.e., as soon as the 
stability argument of fails. The instability is a short distance instability in the sense that 
unstable modes exist for all wavelengths below a certain critical value determined mainly by the 
angular velocity of the hole. However, amongst unstable modes, the shortest wavelength modes are 
the least unstable. 

Having demonstrated the existence of an instability, it is natural to ask what the endpoint of the 
instability is. We propose that a black hole that suffers from this superradiant instability will evolve 
to a stationary, nonaxisymmetric black hole solution. The motivation behind our proposal will be 
explained at the end of this paper. 
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This paper is organized as follows. We start by reviewing the black hole solutions of interest. 
We derive equations governing perturbations of these solutions in section EJ Our strategy for ana- 
lyzing stability is outlined in section 01 The AdS case turns out to be more straightforward than 
the asymptotically flat case, so we analyze stability of AdS black holes in section |21 before analyzing 
asymptotically flat black holes in sectional We discuss the implications of our results in sectional 
We also provide an Appendix, which contains many of the calculational details of the paper. 

Note added: it can be shown that a general gravitational perturbation can be decomposed into 
scalars, vectors and tensors on CP^ , where the tensors satisfy the "doubly transverse" condition 
referred to in the text and the vectors satisfy an analogous condition. Therefore it seems appropri- 
ate to redefine what one means by "tensor" and "vector" to include this condition. The equations 
of motion for scalar, vector and tensor perturbations then decouple from each other without any 
further restrictions. The present paper deals with the most general tensor perturbations. Further 
details will appear elsewhere j21j . 



2 The background solution 

The higher-dimensional generalization of the Kerr solution was obtained by Myers and Perry and 
subsequently generalized to include a cosmological constant in five 8^ and higher ^ dimensions. It 
is parameterized by a mass parameter M and [{D — l)/2] angular momentum parameters Oj. In 
D = 2N + 'i dimensions with equal rotation parameters ai = a the solution is cohomogeneity 1. The 
metric can be written as:^ 

ds'^ = -f{rfdt^ + girfdr"^ + h{rf[di^ + Aadx" - Vl{r)dtf + r^gabdx^dx'' (1) 
,,2 f. 2MH 2Ma^Y^ ,,,2 2 A 2Ma'^\ . 2Ma 



r a?' 
^^''^^ g{r)h{ry ^ = ^"£2' (3) 

where Qah is the Fubini-Study metric on CP^ with Ricci tensor Rah = 2(A^ -|- l)gah, and A = Aadx"- 
is a 1-form such that J = i^dA is the Kahler form on CP^ . This way of writing the metric arises 
from the fact that S"^^^^ can be written as an 5^ fibre over CP^ . The fibre is parameterized by 
the coordinate V') which has period 27r. 

The spacetime metric satisfies R^y = —i^^{D — l)g^i,- Asymptotically, the solution approaches 
anti-de Sitter space with radius of curvature i. An asymptotically fiat Myers-Perry black hole can 
be recovered by taking i oo. The event horizon located at r = r+ (the largest real root of g~'^) is 
a Killing horizon of = dt + Vtud^, where the angular velocity of the horizon is: 

2Ma 
= rf^^ + 2Ma^- 

The mass E and angular momentum J (defined with respect to 9^) are |22j 

where A2N+1 is the area of a unit 2A^ -|- 1 sphere. 

^This can be obtained from the metric of [5] in "unified" Boyer-Lindquist coordinates by changing the radial variable 
r2 l!±fii and M ^ 
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As written the metric is parameterised by (M, a). We shall assume a > 0, which can always 
be achieved by i — >■ — i if necessary. Sometimes it will be convenient to work with more "physical" 
variables (Oif,r+). Fortunately one can easily invert for (M, a) in terms of (riij,r+): 



2{l+rli-^-rlniy l + r2£-2- ^ 

For given r_|_, existence of a regular event horizon imposes an upper bound on iln- 




(7) 

The extremal solution saturating this bound has a regular but degenerate horizon. 

In the asymptotically AdS case, the "co-rotating" Killing vector field is timelike everywhere 
outside the horizon if Qh ^ 1/^ but becomes spacelike in a neighbourhood of infinity otherwise. 
With respect to the metric on the conformal boundary, ^ is timelike if 0,h < l/^j null if = 1/^ 
and spacelike otherwise. For this reason, black holes with Qh > l/i are said to be rotating faster 
than light. Note that the extremal black holes always rotate faster than light, and that such black 
holes can be arbitrarily large compared with the AdS length. 



3 Perturbation equations 
3.1 Scalar perturbations 

Consider a scalar field $ obeying the Klein-Gordon equation 

- = 0. (8) 

We take a separable Ansatz: 

$ = e-*'^*+*™'^i?(r)y(.x), (9) 

where m is an integer. It is convenient to imagine performing a Kaluza-Klein reduction along d/dtp. 
This reduces 5*2^+^ to CP^ . In the reduced description, there is a gauge field A on CP^ with field 
strength 2J, and Y transforms as a scalar field of charge m with respect to this gauge field. We 
assume that Y is an eigenfunction of the charged scalar Laplacian on CP^: 

- p2y = Ay, (10) 

where 

Va = Da- imAa (11) 

is the gauge-covariant derivative on CP^ . The Klein-Gordon equation then reduces to an equation 
for R{r). This can be written in the form of a time-independent Schrodinger equation 

(P 

- -^^{x) + V{x)-^{x) = Q, (12) 

where 

* = ^/hr^R, (13) 

X is a "tortoise" coordinate defined by 
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and the potential V is 



V{x) = Vo-iu;-mnf + f h + !^ + ^] , 



where 



fWhd_ 



fhd_ 
r dr 



{Vhr 



(15) 



(16) 



3.2 Gravitational perturbations 

We consider a hnearized metric perturbation g^i, g^^ + h^^ that is transverse and traceless with 
respect to the background metric: 



flU 



flU 



0. 



The Hnearized Einstein equation is 

AlV = -2^"'(^-1)V> 
where the Lichnerowicz operator is defined by 

= -2R^p,,hP^ - V^h^, - 2l-\D - 1) V- 



(17) 
(18) 

(19) 



Perturbations can be decomposed into scalar, vector and tensor parts according to how they trans- 
form under isometrics of the CP^ base space. In general, these different types of perturbation will 
be coupled together because a perturbation of one type can be converted into one of a different 
type through the operations of differentiation and contraction with the Kahler form. For example, a 
vector quantity Va can be converted into a scalar J'^^[dV)ah- However, we shall identify a particular 
subset of tensor perturbations that can be decoupled from everything else. 
A general separable tensor perturbation takes the form 



0, 



hn 



Re 



(20) 



where Yab{x) is a charge m symmetric tensor on CP^ . Tracelessness of h^y implies that Yah must 
be traceless on CP^ . Transversality of huy is equivalent to 



'D^'Yab 



0, 



(21) 



where P is the gauge-covariant derivative on CP^ defined above. It is well-known that there are no 
tensor modes on CP^ = S'^ so we shall assume N >2 henceforth. 

In Appendix it is shown that the linearized Einstein equation reduces to two equations: 



r'VaYbc = 0, 



and 
= 



f h' 2N\ , , . 

J + - + — ] g'^R' - 



+ 



1 



2im 



J,Y] 



ab 



h r J 

7~Jac{J', Yjcb 



{uj - mflf R + m^h-'R - ^^^^^^R 



R, 



(22) 

Yab 

(23) 



where a prime denotes a derivative with respect to r, and is a gauge-covariant version of the 
Lichnerowicz operator on CP^: 



AtYab = -2RacbdY"' - V'Yab + 4(iV + l)Yab, 



(24) 
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where Rated is the Riemann tensor of CP . 

Equation ((22) arises from the couphng between vectors and tensors: tensor perturbations can 
source vector perturbations and since we have set the latter to zero, the source must vanish, which 
gives equation H22|) . In other words, this equation is the condition that a tensor perturbation must 
satisfy if it is to decouple from vector perturbations. If Yah satisfies both (|21|) and (|22j) then we shall 
refer to it as "doubly transverse". We will demonstrate later that such tensors do indeed exist. 

Define a map J on symmetric two-tensors by 

{JY)ab = Jah^^cd- (25) 

This maps traceless tensors to traceless tensors, and has eigenvalues ±1 whose eigenvectors we refer 
to as hermitian or anti-hermitian respectively. It also maps doubly transverse tensors to doubly 
transverse tensors and commutes with A^. Hence we can simultaneously diagonalize J and on 
the space of traceless doubly transverse tensors. This implies that we can classify eigenfunctions 
of A^ into hermitian and anti-hermitian. Further, in the antihermitian case, one can distinguish 
between Yah of type (2, 0) and (0, 2), which we define by Ja^^A = e^^afe, where e = —1 for (2, 0) 
and (0, 2) harmonics respectively. 

We can now plug all of this into the remaining equation This reduces to an equation for 

i?(r) which can be rewritten in the Schrodinger form (|12|) by defining ^' as in H13() . x(r) by equation 
(|14j) . and the potential is 

= V.-(^- maf + + lii^ + ^ - ^'^ + ^> 7 ^ . (26) 

where u = 1 or — 1 for anti-hermitian and hermitian Yah respectively, A is the eigenvalue of A^, and 
Vb is defined by (|T6|) . 

3.3 Harmonics on CP^ 

We need to know the eigenvalues of charged scalar and tensor harmonics on CP^ . Start with the 
scalars. For = 1 we have CP^ = S"^. Charged scalar harmonics on S"^ have been studied in the 
context of magnetic monopoles |23j . The N = 2 case was studied in |24| . The result for general 
is 

X = l{l + 2N)-m^, l = 2k + \m\, A; = 0,1,2,... (27) 

These charged scalar eigenfunctions on CP^ lift to eigenfunctions of the scalar Laplacian on S'^^'^^ 
where I is the usual "total angular momentum" quantum number.^ 

Recall that there are no tensor harmonics on CP^. Uncharged tensor harmonics on CP^ have 
been considered in j26| I24j . It is straightforward to generalize the analysis of |2^ to include charge: 
see Appendix^ for details. The result is that doubly transverse charged tensor harmonics must be 
anti-hermitian (cr = 1). The eigenvalue spectrum is 

N = 2: X = l{l + A) + 12 - m"^ + Aem, I = 2k + \6 - em\ (28) 



(29) 

em — 2 for em > 4. 

em = 6-l, 6-1 + 2, 6-1 + 4, I, 1 + 2, / = 2,3,4,... (30) 

■^Note that m — mi where (mi, m2, . . .) is the weight vector of S0(2A'' + 2) corresponding to the harmonic in 
question. 
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where 

' 0, 1,2, . . . if em < 4 

k=i 1,2,3... if em = 5 

2,3,4,... if em>6 

This implies that the minimum value of / is /min = 6 — em if em < 3 and /min 
Conversely, for given /, the allowed values of m are given by 



Once again, these harmonics hft to uncharged tensor harmonics on S with / the total angular 
momentum quantum number. Note / > 2, as expected for tensor harmonics. 

Uncharged tensor harmonics on CP^ have been studied in with the result that doubly 
transverse tensors have eigenvalues 

m = 0: X = l{l + 2N)+m + Aa, (31) 

where I = 2k, k = 0, 1,2, . . . is once again the total angular momentum on 5^^"*"^. However, the 
results of should be treated with some caution because this reference overlooks the fact that 
a = —1 is not permissible for CP^, and also gives an incorrect lower bound for /, i.e., / > (whereas 
we know from the above that / > 6 for m = on CP^). The problem with the analysis of [23 is 
that it does not investigate the circumstances under which the would-be eigenfunctions vanish. 

We can obtain a general formula for the charged tensor eigenvalues by comparing our results 
with those for non-rotating black holes. If we set a = then we have a Schwarzschild(-AdS) black 
hole and our tensor perturbations should form a subset of the tensor perturbations considered in 
[110 El- Demanding agreement between our results and those of [Hj yields 

A = /(/ + 2N) +m + Aa-m'^ + 2e(l + a)m. (32) 

This is consistent with all of the above results. The limitation of this approach is that it does not 
tell us which values of / are permitted beyond the obvious restriction / > 2. 



3.4 Unified potential 

To summarize, the equations of motion for both scalar field and (doubly transverse tensor) gravita- 
tional perturbations can be separated and reduced to a single radial equation, which can be written 
as 



+ = 0, (33) 



or, in terms of r, as 

f rl / f Wil/ \ 

+ = 0. (34) 



g dr \g dr 

Using our results for the eigenvalues on CP^ we can present a general expression for the potential 
that applies to both scalar field and gravitational perturbations: 



V = Vo + /V^ -{^- m^f + ^ 



/(/ + 2iV) - ( 1 - ^ ) + 4(1 - a) ( ^ - 1 



(35) 



where Vq is defined by (|l(ij) . For scalar field perturbations, N > 1, a = 1 and I = 2k -\- |m|, 
k = 0, 1, 2, . . .. For gravitational pertubations, N > 2, fi = and the permissible values of cr, / are 
discussed above. 

Note that anti-hermitian (a = 1) gravitational perturbations obey exactly the same equation as 
a massless scalar field. Therefore, one might think that results concerning the stability of (asymp- 
totically fiat) MP black holes against massless scalar field perturbations |181 1281 129j would imply 
stability with respect to a = 1 gravitational perturbations. However, these results concern black 
holes in five jl8j and six [281 129j dimensions with a single non- vanishing angular momentum whereas 
we are interested in black holes in seven or more odd dimensions with all angular momenta equal and 
non- vanishing. It appears that scalar field perturbations of such black holes have not been considered 
previously. Furthermore, there is a difference between scalar field and gravitational perturbations: 
the lower bound on / is different for the two cases. (For N = 2, the lower bound for scalars can be 
either less than, or greater than, the lower bound for doubly transverse tensors, according to the 
value of m.) 
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4 Stability analysis 
4.1 Boundary conditions 

The horizon is located at r = r_|_, which corresponds to x — > — oo: 



x~-logf^^^), (36) 



a \ 
where 

a = — — ^ . (37 

h{r+) 

At the horizon, V — > —{co — rnQH)"^ so the solutions behave as exp (ibi(w — mQf{)x). Regularity on 
the future horizon requires that we choose the lower sign, so we have the boundary condition 

^' = exp {-iiuj - mVtH)x) (38) 

where <l>(r) is smooth at r = r+. Substituting this into the Schrodinger equation and expanding 
around r = r+ gives (choosing <l*(r+) = 1) 

In the asymptotically flat case, we have 

x = r + (r~^^+^) as r ^ oo. (40) 

The general solution as r ^ oo is a superposition of outgoing and incoming waves (proportional to 
gifcr ^-ikr respectively): 



ZontH\%{kr) + Z,^H\%{kr)\ , iV > 2 (41) 



where H^''"'^ are Hankel functions and k = \/ ufi — ■ In the = 1 case, the solution is as above 
except the order of the Hankel functions is now {{I + 1)2 + 2M(/i2 - 2u;2)]i/2. In the asymptotically 
flat case, we are mainly interested in gravitational perturbations for which = and k = uj. 

In the asymptotically AdS case, we have x ~ —(^jr ^ as r — > oo. The asymptotic behaviour 
of the potential is 

^ ~ V^o + ^ ~ ^ ((A^ + 1)^ - ^ + l^H') , (42) 



with corresponding asymptotic solutions ^' ~ r i/2±V(^+^)^+^'^^^. Choosing the normalizable solu- 
tion corresponds to the boundary condition 

^ _ as r ^ oo. (43) 

For stability of the AdS background we demand that a scalar field obeys the Breitenlohner-Freedman 
bound 133] 

>. (44) 

A linearized instability of the black hole would correspond to a solution of the radial equation that 
is regular on the future horizon and vanishing at infinity, with Im(c(;) > 0. In the asymptotically flat 
case this requires Z\^^ = 0. Note that such a solution vanishes exponentially at the horizon. 
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4.2 The case m = 



It is easy to show that there can be no instability (whether asymptotically flat or asymptotically 
AdS) when m = and fJ,^ > 0- If m = then it is natural to consider the potential V = V + lo"^, 
which does not depend on uj. The radial equation becomes 

(P 

- -^^(x) + V{x)^{x) = ^^^(a;) (m = 0). (45) 

Assume that ^ describes an unstable mode, so Im(a;) > and ^ vanishes at the horizon and at 
infinity as described above. The differential operator on the LHS of (|45j) is self-adjoint on such 
functions and hence to'^ must be real so to is pure imaginary and lo"^ is negative. Now consider the 
behaviour of Vq. Letting X = {r /h){d/ dr){^hr^), we have 



r^+i d^ 



It is easy to show that X > and dX/dr > for any > 1. It is also easy to see that [d/dr)g^^ > 
outside the horizon. Hence Vq is positive. If we assume > then the remaining terms in V are 
non-negative so V is positive. Hence o;^ must be positive, which is a contradiction. 

For AdS black holes, it would be interesting to see if this conclusion could be extended to 
tachyonic scalars satisfying the Breitenlohner-Freedman bound by combining our argument with 
that of (33|. 



4.3 Strategy 

We can look for unstable modes using the strategy adopted by Press and Teukolsky for the Kerr 
black hole l3Dj. We expect the black hole to be stable for small angular momentum because we know 
that the higher-dimensional Schwarzschild black hole is stable f^. Hence, for small angular velocity, 
the only admissible solutions of the radial equation must have negative imaginary part, i.e., they 
are quasi-normal modes. If an instability is to appear as the angular velocity increases then one of 
these quasi-normal modes must cross the real axis in the complex a; plane. ^ Hence we expect the 
onset of instability to be indicated by the appearance of a real frequency mode at a critical value of 
the angular velocity. The strategy is to look for such a mode. By continuity it must have Zi^ = 
in the asymptotically flat case, i.e., it must be purely outgoing at infinity. In the AdS case it must 
obey the " normalizable" boundary condition ()43() . 

Note that the radial equation is invariant under lo —u and m —m and, for tensors, e — e.^ 
Hence we can assume lv > 0. Following Press and Teukolsky [321, consider the Wronskian of ^ and 
^ for real to. This must be constant so we obtain 

Im[^'a,^]^^ =0, (47) 

for any xi and X2. Taking xi at the horizon and X2 at infinity and using the boundary condition at 
the horizon and infinity gives, for the (massless) asymptotically flat case, 

(mnH - = ^ (l^outl' - l^inl') . (48) 
Hence a purely outgoing mode must have 

< a; < m^H (asymptotically flat). (49) 

^Various mathematical subtleties such as modes coming in from infinity might invalidate this statement but such 
subtleties do not occur for Kerr |31| and we shall ignore this possibility here. 

^This just corresponds to complex conjugation of the solution. This invariance arises from a discrete symmetry of 
the background which, in the coordinates of |9j is f — > —t, (pi —(pi. In our coordinates this amounts to t —> —t, 
■0 —ip and A ~A. 
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In other words, the mode must be superradiant. In the AdS case, the LHS is unchanged but the 
term at infinity vanishes and we must have 

uo = mQ.H (asymptotically AdS). (50) 

Physically, this is simply the statement that there cannot be a constant flux of radiation through the 
horizon if the flux at infinity vanishes. Note that for both cases, we must have m > since > 0. 

We shall discuss the asymptotically AdS case first because the fact that we only have to consider 
a single value of uj makes this case simpler to analyse than the asymptotically flat case, for which 
we have to consider a range of values for u. 

5 Asymptotically anti-de Sitter black holes 

5.1 Form of the potential and behaviour of solutions 

Set UJ = itiUh > 0. The potential vanishes at the horizon and is monotonically increasing just 
outside the horizon. For large r, the potential increases^ proportional to r^. What happens in 
between depends on the values of the parameters. For small uj, V is positive everywhere. However, 
for sufficiently large to, there is a "classically allowed"^ region in which V is negative. In more detail: 

V has roots at r = ri , r2 . The potential is positive for r+ < r < ri , negative for ri < r < r2 and 
positive for r > r2, i.e., there is a potential barrier separating the classically allowed region where 

V is negative from the horizon. 

Note that, in the AdS case considered here, the initial data ^'(r_|_) and \I''(r_)_) are real and 
positive (since u = itiUh)- Hence ^ is real everywhere. It is easy to see that ^ will simply increase 
monotonically if V is positive everywhere. Hence we need V to be negative somewhere for an 
acceptable solution of the radial equation to exist. The solution ^ will increase monotonically in the 
potential barrier, oscillate in the classically allowed region, and then match onto a sum of growing 
and decaying^ terms at large r. We need to tune until the coefficient of the growing mode 
vanishes, i.e., until we obtain a "bound state" solution of the radial equation. 

5.2 Small AdS black holes 

Consider the case <C i. For r ^> r-|_, the presence of the black hole does not affect the potential, 
which is determined entirely by the AdS background. For r ^ i the potential increases as and 
for r+ r <^ £, there is a centrifugal barrier in which F ~ This matches onto the decreasing 

section of the potential following its local maximum at r ~ r+. If a; = mQu is large enough then V 
will be negative at its local minimum, which occurs at r ~ 

The solution in the region r 3> r+ is easy to obtain since we can approximate it by a solution in 
pure AdS. The solution obeying " normalizable" boundary conditions as r — oo is 

* = -'^^^^/^ (l + ^ j F[a,c-b;a-b+l; , (51) 

where F is the hypcrgcomctric function, and 

uoi + l + N + 1 + ^{N + 1)^ + fi'H'^ ^_ Loe + l + N + l- V(jV + 1)^ + 

" ~ 2 ' ^" 2 ' 

^If /u^ is close to the Breitenlonher-Preedman bound then the coefficient of proportionaUty is negative but we shan't 
worry about this and our results for small black holes suggest that it doesn't change the qualitative behaviour of 
solutions. 

^Of course, everything we arc doing is classical but since we have written the radial equation in the form of a 
Schrodinger equation, we can borrow terminology such as "classically allowed", "bound state" etc. from quantum 
mechanics. 

'^More precisely: non-normalizable and normalizable. 
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c = 



l + N + l. 



(52) 



Now pick ro so that r+ ^ tq ^ i. Near r = rg, the solution (|5H1 is a superposition of a growing solu- 
tion ~ 7-^+^+1/2 a, decaying solution ~ j.i/'^-i-N ^ T\ie latter will dominate for sufficiently small 
ro (ro can be made as small as we like by decreasing r+). Hence the solution will be monotonically 
decreasing near r = ro. However, r = ro is inside the potential barrier outside the horizon of the 
black hole, where we know that ^ should be increasing! So in fact the coefficient of the decreasing 
solution must vanish, which leads to the quantization rule c — b = —p where p is a non- negative 
integer. This implies lo£ = 1 + N + 1 + y/ {N + 1)-^ + ^^^2 _|_ 2p. This is just the usual quantization of 
frequency in global AdS. However, here we know that uj = rnQn- Hence, for small AdS black holes, 
a real frequency solution to the radial equation exists for 

a„i^:±ii±i±^m±wT^l±^^ p=o,i,2,... (53) 

m 

The solutions with p > correspond to "excited states" for which the solution of the radial equation 
oscillates (with p + 1 extrema) before approaching zero at large r. We are interested in the onset of 
instability, corresponding to the smallest value of ^Ih^ for which a solution exists, so we are mainly 
interested in p = 0. 

For scalars, we have I > m so, for given m, the smallest value of ^Ih ^or which we have a solution 

is 

m 

This is the critical value of 0,h beyond which modes with angular quantum number m become 
unstable. Note that it always exceeds consistent with the proof of stability for < 1/i 
given in JHI- However, this proof has been criticized j34j because it assumes the dominant energy 
condition, which is violated if < 0. Our result shows that, for small black holes, this does not 
matter so long as the Breitenlohner-Freedman bound is satisfied. 

Note that the critical value for Qh tends to 1/i from above as m ^ oo. This proves that, 
for small black holes at least, the instability sets in as soon as exceeds l/l, with the shortest 
wavelength modes becoming unstable first. 

Now consider gravitational pertubations, for which the threshold of stability occurs at 

n„, = i±^. (55) 

m 

For N = 2, taking / = Zmin("T') and e = 1, this evaluates to 11,5,3 for m = 1,2,3 and 1 + 4/m for 
m > 4. (Taking e = — 1 just makes Imm bigger.) So the conclusion is the same as for scalar field 
perturbations: a superradiant gravitational instability sets in as soon as the angular velocity exceeds 
the speed of light, with the shortest wavelength modes becoming unstable first. 

In the context of the AdS/CFT correspondence, we have to consider a higher dimensional space- 
time consisting of the product of the black hole space-time with a compact internal space (e.g. a 
sphere). It is believed that small AdS black holes are unstable with respect to the Gregory-Laflamme 
instability under which they are expected to localize on the internal space. This means that, 
although small black holes with 0,Hi < 1 do not suffer from a superradiant instability, they are 
nevertheless unstable. In order to eliminate the GL instability we have to extend our results to large 
AdS black holes, i.e., r+ > i. This will be addressed in the next section. 

5.3 Numerical analysis 

Fix r_(_ and {l,m). For given we integrate numerically the equation for ^ (with uj = mQH) out 
from r = r_|_. For small i^n^ we find ^I' ^ cxd as r — > oo. However, when a critical value of 0,Hi is 
exceeded, we find ^ —oo as r ^ cxd. The critical value corresponds to a solution for which ^' ^ 
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Figure 1: Plot of solution \l'(r) against r jl for doubly transverse tensor perturbations with N = 2, 
l = Q^m = 8, r+/i = 0.7, ^ni = 1.46242. 

as r — > oo, which is what we are looking for. The precise value of ^Ih^- at which this occurs can be 
obtained accurately by interval bisection. This can then be repeated for different values of r^.., I, m. 

In this section, we shall consider only gravitational perturbations with N = 2 (D = 7). A typical 
numerical solution for the wavefunction at a critical value of ^Ih is shown in figure ^ Results for 
the lowest values of m are shown in figure |2 We plot the critical value of Qh as a function of r+ for 
different values of m with / = /min("i-) and e = 1. These curves terminate when they hit the curve 
corresponding to the extremal black hole (beyond which no black hole solution exists). There are 
two points to note about these results. First, as r^/i — > 0, they are in excellent agreement with the 
theoretical predictions for small black holes discussed above. Second, the only effect of increasing r+ 
is a very slight decrease in Qh- This decrease was too small for our numerics to detect for m = 1,2. 
For (/, m) = (3, 3), the critical value of ^Ih decreases from 3.00 to 2.99 as r_|_ increases from to its 
maximum. For {l,m) = (2,4) the decrease is from 2.00 to 1.85 and for {l,m) = (6,8) it is from 1.50 
to 1.42. 

To obtain results for r^/i > 1 we have to go to higher m. Some sample results are presented 
in figure 01 Our numerical results suggests that the critical value of ^Ih^ is always greater than 1, 
and tends to 1 as m ^ oo, i.e., the curves on the figures accumulate at i^n^ — *■ 1+ as m ^ oo. 
The interpretation is exactly the same as for small black holes: for any given r+, if we start from 
rin = and increase i},H then as soon as i^n^ exceeds 1, the black hole will become unstable to all 
perturbations for which m exceeds some critical value. In the next section we shall demonstrate this 
analytically. 

5.4 WKB analysis 

In this section we consider both scalar field and gravitational perturbations governed by the effective 
Schrodinger equation with potential (|35|) . As before, we are interested in modes at the threshold of 
instability so a; = mQu- The strategy is to look at the potential for large m (and hence large lo). 
In this limit, we can use WKB techniques. 

As discussed above, if an appropriate solution of the radial equation is to exist then the qualitative 
behaviour of the potential must be as follows: it vanishes at the horizon r = r_(_, increases to a local 
maximum and then decreases and becomes negative in some range ri < r < r2. For r > r2 the 
potential is positive and monotonically increasing. 

We are interested in bound states, for which ^ obeys the boundary condition (|43() . The appro- 
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Figure 2: Plot of Q^ni against r+/£ for doubly transverse tensor perturbations with N = 2, 
e = 1. No black holes exist in the empty region bounded by the curve in the top right 
of the diagram. Curves rij^(r+,i,m) corresponding to existence of a normal mode with fre- 
quency u = mClH are displayed for different values of l,m. Prom top to bottom, {l,m) = 
(5, 1), (4, 2), (3, 3), (2, 4), (3, 5), (4, 6), (5, 7), (6, 8). 
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Pigure 3: Plots of critical value of Qh for doubly transverse tensor perturbations for iV = 2, e = 1, 
{l,m) = (38,40) (top) and (78,80). No black holes exist in the empty region bounded by the curve 
in the top right of the diagram. The fact that the Qh curves do not quite meet this curve is due to 
the limitations of our numerical method. Note the scale of the vertical axis. 
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priate WKB solution in the region r > r2 is 



X2 



^(x) « 7"^/2 exp ( - / jdx ) (56) 



where 7(x) = y/V{x). Note that, for large r, we have 



^2 

V{x)^^, x^O- (57) 



where 



K = y(A^ + l)2-i + /i2^2. (58) 

The WKB approximation is valid when \^~'^d'y/dx\ and \^~''^d'^'y /dx"^] are small, which requires 
K ^ 1. Hence we do not expect our results to be reliable for a scalar field with mass close to the 
Breitenlohner-Preedman bound. For gravitational perturbations fi = so a priori we might not 
expect WKB to work very well. However, the same remark applies to the calculation of black hole 
quasi-normal modes, where WKB has been found to be accurate [23 • We shall see that the WKB 
results are in good agreement with our numerical results for ^ = 0. In any case, the WKB method 
is certainly reliable for sufficiently massive scalar fields. 

The standard connection formula (see, e.g. jSH]) imply that, in the classically allowed region 
ri < r < r2, we have 

^ « 2k-^/^ cos (^y""' kdx - ^) = 2/c-^/2 J^j _ J"" _ ^ _ (59) 
Here k{x) = \/—V{x) and 

1=1^ kdx. (60) 

J Xl 

Applying the connection formulae once again leads to the WKB solution for r < ri: 

^' 27"^/^ cos /exp ( / -ydx - f -fdxj . (61) 

\J—oo J ~oo J 

The WKB approximation breaks down as x ^ —00 (r — > r-|_) because ^ — > 0. This requires separate 
analysis. Using and (|H7|) and setting y = exp(ax/2), the Schrodinger equation becomes 



dy^ dy \ 



y -FT + y-r-[ 2 y + ^(r) pi' = o- (62) 



This is Bessel's equation for y — > 0. The solutions are Iq{z) and Kq{z) where 



z=-Jr+V'{r+)y. (63) 

rv V 



As m is large, V' is large and hence for moderate y, the argument of the Bessel functions is large 
and we can approximate them by their asymptotic forms: 

Uz) ~ K,{z)^ J^e-\ (64) 

In general we have a linear combination of these two solutions. However, if ^ is asymptotically 
decaying for large z then it must be proportional to Kq. This is precisely what we have found from 
our WKB analysis: equation (|6H) decays for large z. Hence our WKB solution matches onto the Kq 
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solution near the horizon. However, the Kq solution diverges at the horizon, and the only way to 
avoid this is if cos / = 0, i.e., 



p+-j7r, p = 0,l,2,... (65) 
The next step is to evaluate 

1= r J-V{r)Ur. (66) 

We assume (motivated by our numerical results and results for small black holes) that = 
1 + 0{l/m) and that the roots ri and r2 are of order y/ml. These assumptions can be justified a 
posteriori. We take m ^ 1, m ^ r^/^^. Then within the integral r is of order ^/m^ and we can 
approximate g/ f ^ ^^/r^. Now consider H35|l with the above restrictions. We find 

v{r) . p + ™^ (r- - al + ,■-) + + (67) 

We assume I ~ /min = m + 0(1) for large m. Note that we must have > 1 for F to be negative 
somewhere. The roots ri and r2 are determined by V{ri) = 0, with solutions 



We then find 



2 ' - l) - 2m(/ -m + N)±^{m? {p.]jP - l) - 2m{l - m + N)Y - Ak^t 



m" \l„r - i - imu - m + iy)± \ [m^ - i - ImiL - m + iV - 4K^m 



(68) 

^ / " ^ V(r2 - rl){rl - r^) = ^ ' , (69) 
Jri 4rir2 

where the integral was calculated by changing variables to = and working in the complex 
ty-plane.^ Plugging this into the quantization condition H65j) gives, at large m with I — m = 0(1), 



^Hi-l + ^\ l-m + N + 1 + ^{N + IY + /i2^2 + 2pj . (70) 

Finally we can use these values for Q^h to verify that rj are really of order y/ml, as previously 
assumed. Recall that validity of the WKB approximation requires that {N + 1)^ + //^^s ^ -^j^ 
this result agrees with our result (|53|) for small black holes. Note that our result is independent of 
r_(_, which explains (at least for large m) why our numerical results showed little variation with r_(_. 
Setting /u = 0, = 2 and p = 0, the WKB method gives ~ 1 + 3.96/m in good agreement with 
our numerical results so, as advertised above, the WKB method works well even though one might 
question its validity when ^ = 0. 

The coefficient of 1/m in (|70)) is positive (at least in all cases for which we know /min) so once again 
we see — > 1+ as m — > c« so, irrespective of the size of the black hole, once its angular velocity 
exceeds the speed of light it becomes unstable to perturbations of arbitrarily short wavelength. 



5.5 WKB calculation of unstable modes 

So far, we have been looking for real frequency modes, whose existence indicates the onset of insta- 
bility. In this section, we shall look for unstable modes directly using the WKB approach. This has 
been used previously in a study of the superradiant instability of the Kerr black hole in the presence 
of a massive scalar field |12j . 

We allow oj to be complex: oj = ujr + iujj and assume \ujj\ <C = 0(m). Neglecting a;/, the 
potential is now negative at r = r+ and has an additional zero at r = ro, ?"+ < tq < ri. The 

*More precisely, introduce a branch cut between rf and r^; the integral around the branch cut can be deformed 
and gets contributions from a double pole at w = and the circle at infinity. 
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potential is negative in a region r_|- < r < ro near the horizon, positive in the potential barrier region 
ro < r < ri, negative for ri < r < r2 and positive for r > r2- It is straightforward to extend the 
WKB analysis to cover this potential so we shall neglect the details. Assuming ujr < inQn and 
I — m = 0(1) we obtain the quantization condition 



+ l (/_m + iV + l + A/(A^ + 1)2-1 + ^2^2 + 2p) , p = 0,1,2,... (71) 
m m \ V 4 I 

The WKB solution for r > ri takes the same form as before. For vq < r < ri the solution is (7 = 
again) 

^ « (-)^7"^/^ exp J^' -fdx^ exp -fdx^ . (72) 

For r+ < r < ro the solution is {k = ^/—V) 

^ PS {-fk-^/^ exp -fdx^ exp - i £° kdx^ . (73) 

To calculate uji, write the Schrodinger equation as 

- dl"^ -{lu- mVtf ^ + = 0, (74) 

where V is independent of lo. Multiply this equation by ^ , take the imaginary part and integrate 
(wrt x) from r = r,, > r_|_ to r = cxd to obtain 

rr=co _ 

2uJi dx{ujR-mn)\^\^ = lm[^d^^]^^^^. (75) 

Take < ro and evaluate the RHS using the WKB solution. The integral on the LHS is dominated 
by the classically allowed region ri < r < r2 where ujji ^ mfl. We obtain 

" 2.^/-|^|2d. • (^6) 

This is the same as equation 41 of 12 . The denominator can be calculated using the WKB solution 
and the expression for ujr. It turns out to be 0(1) at large m. In the numerator, there is a 
contribution of 0{m) to the integral coming from the region near r = ro so we can bound 

< ujjl < aex.p{-f3m) (77) 

for some positive constants a, (3. This justifies our assumption a;/ <C ujr. We see that, although large 
m modes are the first to become unstable when VLh^ exceeds 1, the growth time of the instability is 
exponentially large in m so these modes are the least unstable. This suggests that the most unstable 
modes will be those for which m is not particularly large. It would be interesting to calculate a;/ for 
such modes. 



6 Asymptotically flat black holes 
6.1 Introduction 

We know that a mode at the threshold of instability must obey < a; < rnQn, i-e., it is superradiant. 
The only known way that superradiant modes can actually lead to an instability is if they can be 
trapped by the potential at infinity, i.e., they must be bound states. This would require a local 
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minimum in V, as in the AdS case, or for a massive scalar field in four dimensions. In all cases 
that we have examined, the qualitative form of the potential for /x = is: V ^ —{lo — ni^H)'^ as 
X — > — oo, then V increases to a positive maximum and decreases to — u;^ as x ^ oo. It appears 
that a local minimum in V is not possible so there is no obvious sign of any gravitational instability 
apparent from our radial equation. 

This qualitative argument is no substitute for a quantitative study. We shall analyse the radial 
equation in two cases: first for large m and lv using the WKB method and then numerically for 
D = 7. 



6.2 WKB approximation 

Consider large m with r+w/m fixed and / ~ /min so l/m — > 1. We find 

W= lim ^ = - ^\ (78) 

This is true even with a mass term (as we're considering a limit in which w ^ /x). The explicit form 
of is 

Pir) 

W = ('79) 

where 

p(^) = _i^^2iV+2 + ^2iV _ fl - ^ V . (80) 

The qualitative behaviour of P[r) is increasing from a negative value at r = r_|_ to a maximum at 
= Nm? /{[N + l)i^^), then decreasing to — oo as r — > co. The important question is whether the 
maximum of P is positive or negative. If it were negative then V would be everywhere negative 
so the WKB solution proportional to exp(i / ^J—V dx) would be valid everywhere. This solution 
would be regular on the future horizon (for to < mQ.H) and purely outgoing at infinity, i.e., the 
WKB approximation would give Z^^ = 0. We would then have to examine whether this conclusion 
persisted beyond the WKB approximation. If so, we would have evidence of an instability. However, 
in Appendix 10 we show that this does not happen: the maximum of P is positive in the range of 
interest < a; < ttiVLh- 

Qualitatively, the form of V is as follows. It takes the value —{uj — milH)'^ at the horizon, 
increases to a positive maximum and then decreases to — a;^ at infinity. In other words, there is a 
potential barrier of height proportional to m? separating the classically allowed region near infinity 
from the classically allowed region near the horizon. The WKB method will then give |^out/-^in| ^ 1* 
However, as argued above, a mode at the threshold of instability will have |Zout/-^in| oo. We 
conclude that no such mode exists for large m and to. 

The degenerate case in which P vanishes at its maximum occurs only for extremal black holes 
with u) = mil.H- Equation ()48|) then gives |Zout/'^in| = 1, so there is no evidence of any instability 
in this case either. 



6.3 Numerical results: asymptotically flat case 

We shall only consider gravitational perturbations in D = 7 so ^ = 0, = 2. 

For given l,m, our strategy (following jHO]) is to start with small and search the interval 
< UJ < mfln for a solution of the radial equation that is regular on the future horizon and outgoing 
at infinity. This is then repeated for increasing values of ^Ih up to the maximum value. 

®In higher dimensions, it appears that even a mass term for a scalar is not enough to lead to a superradiant 
instability, at least for MP black holes with a single non- vanishing angular momentum |39| . 
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To perform the numerical integration of the radial equation, we start by using the variable ^> 
defined by equation (|38() . which satisfies the equation 



- fg-'-^ (f9~'T) + - ^^H)fg-'^ +\v + {u;- mnn?] ^ = 0. (81) 
ar \ ar / ar ^ J 

The advantage of working with ^ rather than ^ is that <I> is a smooth function of r at r = r+. 
Equation H39|) determines the values of ^ and at r = r+. The equation for $ is integrated 
numerically out to r = 2r+, where we then calculate the corresponding values of ^ and }^ We 
then take these as initial data for the integration of the equation for ^ out to large r. 
A convenient object to consider is the ratio 

Z.^=limi^^±^. (82) 

We are looking for a; for which this ratio diverges, corresponding to a purely outgoing solution. We 
evaluate the right hand side at r = pr_|_, p = 10, 11, . . .. Denote these values by Z^. The numerical 
integration is terminated when the sequence Zp has converged to an accuracy of 0.1%. The solution 
is rejected if the Wronskian \mi^dx^^ varies by more than 0.1% during the numerical integration. 

In all cases we have examined, the qualitative form of the potential is the same as we found in 
the WKB analysis above, i.e., a potential barrier with a positive maximum separates the classically 
allowed regions near the horizon and far from the black hole. The corresponding behaviour of the 
solution ^ is: oscillation near the horizon, exponential growth in the potential barrier region and 
then oscillation out to infinity. If the potential barrier is large then this implies that the amplitude 
of oscillation far from the black hole will be large. However, we are looking for a mode with Z\^ = 0. 
From such a mode obeys 

hm = . (83) 

Hence, unless w is very small, such a mode will not have a large amplitude. Hence it seems very 
unlikely that we will find a suitable mode when the potential barrier is large. Phrasing the argument 
slightly differently, if ^ is large then and/or Z^^t must be large compared with the left hand 
side of equation (|48() . which implies Z ~ 1. More physically, if the potential barrier is large then one 
expects almost perfect reflection and very little transmission, so the amplitude of ^ is much greater 
far from the black hole than near the horizon. 

This argument suggests that we should examine the case for which the potential barrier is 
smallest. The potential barrier is minimized when I is as small as possible and m as large as 
possible. The most favourable case (using ((201)) is therefore likely to be / = 2, e = 1, m = 4. 
Our numerical results are shown in figure 0J where we plot Z against uj/{mO,H) for i^H/^max = 
0.5,0.7,0.9,0.99,0.999 where r^max = \/N/{N + l)/r+ is the upper bound on Qh- The curves have 
the same qualitative shape as for the Kerr black hole jHOj, i.e., Z is very close to 1 for small w, then 
increases to a maximum near uj = mfln and decreases back to 1 at a; = m^n (the latter property 
follows from equation H48|)). The position of the maximum tends towards uj — Tfi^u as — ^ ^max* 
The largest value for Z \s Z = 1.115 so there is no sign of Z diverging anywhere, as would be 
required for an instability. Note that the amplification of energy flux in superradiant scattering is 
given by Z^ so the maximum amplification apparent in our data is about 24%, and is achieved as 
the black hole tends to extremality and lo m^n—- This is just as for Kerr, although for Kerr, 
the maximum amplification is much greater: 138% |32j . 

We have repeated our analysis for other values of (/,m). The results are qualitatively similar to 
the case we have just discussed. For ^H/^raa.^ = 0.99, the largest value of Z obtained for e = 1 
and (/,m) = (3, 5), (3, 3) was 1.056,1.000 respectively, reflecting the fact that decreasing m tends 

^"^In practice, we use the freedom to rescale ^ by a phase to set = $ at r = 2r+. 
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to increase the potential barrier. In figure El we exhibit how Z varies with / with m = mniax(0 — 
I + 2. The largest value of Z occurs for the / = 2, m = 4 case discussed above, and Z decreases 
monotonically to 1 as / increases, in agreement with our WKB analysis. 




0.0 0.25 0.5 0.75 1.0 

Lo/{mflH) 



Figure 4: Plots of Z againt u;/(mr2j/) for (from bottom to top) nn/^ms^x = 0.5,0.7,0.9,0.99,0.999 
with e = l,/ = 2,m = 4. Note that Z = 1 for uj = mfln but this point has been deleted from the 
topmost two curves to make the figure clearer. 
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Figure 5: Plot of Z against Z = 2, 3, . . . for m = mmax = / -|- 2, e = 1, QH/^ma.x = 0.99, uj/{mVLH) = 
0.99. 
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7 Discussion 



We have shown that there exists a class of higher-dimensional rotating black hole solutions for which 
linearized gravitational perturbations can be studied analytically. We have concentrated on the 
particular case of doubly transverse tensor modes. The equations of motion for such modes decouple 
from other types of perturbation and, after separation of variables, reduce to a single ODE governing 
the behaviour in the radial direction. It will be interesting to see whether the equations of motion 
governing other types of gravitational perturbation can also be reduced to ODEs. 

Our radial equation could be used to calculate quasi-normal modes of these black holes. This 
might be simpler than the calculation of quasi-normal modes of a Kerr black hole because for Kerr, 
the eigenvalues of the angular equation have to be determined numerically and then fed into the 
radial equation whereas for cohomogeneity-1 black holes, these eigenvalues (our A) can be determined 
analytically. 

In the asymptotically flat case, our WKB argument reveals that, for large m, the gravitational 
perturbations are stable. We have also investigated the D = 7 {N = 2) case numerically for smaller 
TO, finding no sign of any instability. (A similar treatment for D > 7, will need to answer two 
questions regarding charged doubly transverse tensor harmonics on CP^ with N > 2, namely what 
is the lower bound on / and do hermitian (cr = —1) modes exist?) It would be interesting to repeat 
our analysis using the continued fraction technique [HE] as this may prove more efficient. It is also 
of obvious interest to investigate whether the analytic proof of mode stability for Kerr could be 
adapted to our radial equation. 

Our proof that to = modes are stable is relevant to the conjectured existence j41j of higher- 
dimensional rotating asymptotically flat black holes with fewer symmetries than any known solution. 
The motivation behind this conjecture is that a stationary black hole must be axisymmetric (which 
has now been proven for higher dimensional black holes 02), i.e., it must admit a rotational symme- 
try. However all known higher-dimensional black hole solutions have multiple rotational symmetries 
and therefore more symmetry than one would expect of a general stationary black hole. So perhaps 
there exist undiscovered solutions that, although axisymmetric, possess fewer symmetries than any 
known solution. If such solutions are continuously connected to the known solutions then one might 
expect to see evidence for them in perturbation theory as stationary (iv = 0) axisymmetric (to = 0) 
zero modes separating stable low angular momentum black holes from unstable high angular mo- 
mentum ones. However, we have seen that, for doubly transverse tensor perturbations, to = modes 
are stable and there are no stationary axisymmetric perturbations. Hence there is no evidence for 
the conjecture of [U in the sector of doubly transverse tensor perturbations. It will be interesting 
to see what happens for other types of perturbation. 

We have studied the case of asymptotically anti-de Sitter black holes in some detail. Our main 
result is that the superradiant instability of such black holes discussed in JJJI occurs precisely when 
the angular velocity of the black hole exceeds the speed of light (in the sense that the co-rotating 
Killing field becomes space- like on the conformal boundary). In other words, the instability occurs 
precisely when the stability argument of ^U] fails. Our results for small AdS black holes also enable 
us to address a loop-hole in the proof of TO^, namely that it could be invalidated by the presence of 
a tachyonic scalar field (obeying the Breitenlohner-Freedman bound), which violates the dominant 
energy condition. For small black holes at least, such scalars do not behave any differently from 
more massive scalars, i.e., the threshold of instability is unaffected. It would be interesting to see 
whether the methods of jlUj and '35' could be combined to give a more general proof of this. 

Something that has not been addressed in the literature is the end-point of the instability of 
For the super-radiant instability of a Kerr black hole in the presence of a massive scalar, the evolution 
is clear: the black hole gradually loses energy and angular momentum to bound superradiant modes 
of the scalar field. These couple non-linearly to gravitational waves (and thereby to non-super radiant 
modes of the scalar field) so this energy and angular momentum is ultimately radiated to infinity. 
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After a very long time the black hole will have lost all of its angular momentum this way. 

In the AdS case, the evolution of the instability must be rather different. A black hole that is 
initially rotating faster than light will lose energy and angular momentum into superradiant modes 
of all fields of the theory under consideration. However, in AdS this cannot be radiated to infinity so 
instead the energy and angular momentum in fields outside the black hole must accumulate over time 
and backreaction will become important. If the system ultimately settles down to an equilibrium 
state then this must be described by a new stationary, asymptotically AdS black hole solution. 
Presumably the angular velocity of this new solution will not exceed the speed of light. Note that 
this argument does not depend on the details of the theory: it would be valid if gravity were the 
only field present, so there should even exist new vacuum black hole solutions. 

Further evidence for the existence of new solutions comes from our analysis of modes at the 
threshold of instability. Since these have oo = rnQn, they are preserved by the co-rotating Killing 
vector field ^. However they break the symmetries generated by d/dt and d/dip. The existence of 
these modes could indicate the presence of new branches of solutions that bifurcate from the solutions 
of IHlini at the critical (m-dependent) value of ^h- The new solutions would not be invariant with 
respect to d/dt or d/dip. In other words, they would not be stationary or axisymmetric. However, 
^ would still describe a symmetry. Since the bifurcation point occurs when the original black hole 
is rotating faster than light, ^ would be spacelike near infinity but timelike near the horizon. So, 
near the bifurcation point, the new branch of solutions would correspond to black holes that are not 
stationary in the usual sense but nevertheless admit a Killing field that is timelike near the horizon. 
This Killing field becomes spacelike outside an "ergosphere" (this is what happens for the co-rotating 
Killing field of the Kerr black hole). However, if such solutions do exist, they are themselves rotating 
faster than light and therefore seem very likely to be unstable. 

What happens as one moves further away from a bifurcation point? Obviously we can only 
speculate, but one possibility is that, if one moves sufficiently far along one of the new branches of 
solutions, one reaches solutions for which ^ is timelike everywhere outside the horizon. These would 
correspond to genuinely stationary black holes which are nevertheless nonaxisymmetric. There would 
be no violation of the theorem that a stationary black hole must be axisymmetric 42 because this 
theorem assumes that the stationary Killing field is not normal to the event horizon whereas ^ is 
normal to the horizon of all the black holes we have been discussing. If such black holes exist then it 
is natural to guess that these should be the new solutions describing the endpoint of the superradiant 
instability. 

The possibility of a black hole being stationary with respect to a Killing field that does not 
approach the "usual" generator of global AdS time translations deserves further comment. Consider 
AdS^ in global coordinates: 

ds^ = - ^1 + ^ j dt'^ + dr^ + [dO^ + sin^ ed(j)j + cos^ 9d^fj . (84) 

The generator of global time translations d/dt is a globally timelike Killing field whose norm diverges 
at infinity. Now consider the Killing field 

d_ l_d_ 1^ 

dt^ id(Pi^ id(P2' ^ ^ 

This is globally timelike with constant norm. The same construction works in any odd-dimensional 
AdS spacetime. Therefore there are (at least) two qualitatively different ways that an odd-dimensional 

^^This is reminiscent of what happens for black strings. These suffer from the Gregory-Laflamme instability |35| . 
They admit a static zero mode ^43; that signifies the onset of instability |44| and suggests the existence of a new family 
of solutions bifurcating from the known ones . This is indeed the case |4fcil I47| but the new solutions have too low 
an entropy to be the endpoint of the original instability, and seem likely to be unstable themselves. 
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asymptotically AdS space-time can be stationary: the generator of time-translations could have ei- 
ther unbounded norm or bounded norm. This does not appear to have been discussed before. The 
known AdS black hole solutions are stationary in both senses because they admit Killing fields that 
have the same asymptotic behaviour as d/dt and d/d(j)i above. However there may well exist AdS 
black holes with less symmetry that are stationary only with respect to a Killing field of bounded 
norm. 

We do not know whether the stationary nonaxisymmetric black holes discussed above must be 
of this form. If they are, then, since the stationary Killing field must be normal to the horizon, and 
since a Killing field of bounded norm is null on the conformal boundary, such black holes must be 
rotating at the speed of light. Superficially, this makes sense because the superradiant instability 
"switches off' when the rotation of a black hole no longer exceeds the speed of light so one might 
expect the black hole to evolve to a final state rotating precisely at the speed of light. However, the 
evolution of the instability is a time-dependent process, during which the notion of angular velocity 
is not even defined, so we see no reason why the final time-independent state could not be rotating 
slower than light. 

Finally, we note that super symmetric black holes exist in AdS^ |3H1- Supersymmetry guarantees 
the existence of a non-spacelike Killing vector field that is normal to the event horizon [ITj and 
timelike with bounded norm at infinity [IHI) i-e., these solutions have precisely the behaviour that 
we have just discussed and hence rotate at the speed of light. However, these solutions admit extra 
Killing fields analagous to d/d(j)i (and hence also d/dt) above so they are also stationary in the 
usual sense (i.e. with respect to a Killing field of unbounded norm). The existence of these extra 
Killing fields appears unrelated to supersymmetry, which raises the question of whether there exist 
more general supersymmetric black hole solutions without these extra symmetries. Such black holes 
would be nonaxisymmetric, and stationary only in the new sense that we have been discussing. 
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A The Lichnerowicz equation 

Introduce an orthonormal basis 

= fdt, e^=gdr, e'^ = h [dtp + A - Qdt) , e' = r&Jx'', (86) 

where e* is a vielbein on the base. The dual basis defined by (e", ep) = 5p is: 

eo = {dt + ^d^) , ei = g'^^dr, 62 = h~^d^, ei = r'^iti - (A, ei)d^) (87) 

where Cj are the dual vectors to &. The spin connection is defined by (a,/3 are tangent space, [i^v 
are curved indices) 

and obeys 

de" = -(j"^ A e^. (89) 

^^The terminology here may be a little confusing. If one defines angular velocity for such black holes in the usual 
way then it will vanish because the stationary Killing field is normal to the horizon. However, this is angular velocity 
defined with respect to a stationary bulk observer. Since the stationary Killing field is null on the conformal boundary, 
such an observer actually rotates at the speed of light with respect to the boundary and hence so does the black hole. 
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We find 
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where i,j denote components with respect to the base space vielbein, and uj^ij are the components 
of the base space connection with respect to the base space vielbein. 



The curvature two forms are given by i?"^ = duj^'g + uj"^ A ui'^g. We find 
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The only non-zero components of A^/i are (using ()2ipl 
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Plugging this into the 2i component of the linearised Einstein equations leads to equation (|22() and 
the ij component gives ()23p . 
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B Charged tensor harmonics on CP^ 



To calculate the spectrum of A;^ acting on the space of charged doubly transverse traceless tensors 
on CP^ we will use the method of We will use the notation and conventions of j26j (with e = 1) 
for ease of comparison and then convert to our notation at the end. We consider a compact Einstein 
Kahler manifold of complex dimension two, with Kahler metric Qab which satisfies Rab = ^dab and 
has Kahler form J. A gauge field A is defined by dA = (A/ 2) J. 

In it is shown that a general symmetric two-tensor can be decomposed as: 

Ya, = JMb + L^Hb + K'^cb (93) 

where u!^ and uj^ are uniquely determined anti self-dual two forms (the Kahler form J is taken to be 
self-dual), and L is a particular (0, 2) complex self-dual two-form and K = L. Explicitly, in terms of 
complex vierbeins z"^ (so g = z"^ z"^), L = —iz} A z^ and J = ^{z^ A Z"*^ + A z^). The covariant 
exterior derivative on a p-form u) of charge n (with respect to A) \s D = d — inA where d is the 
ordinary exterior derivative . The two- form J is uncharged, L has charge n = —2, K has charge 
n = 2, and all three are gauge covariantly constant (i.e. annihilated by P = V — inA). One can split 
the operator into (1, 0) and (0, 1) parts as D = + where = d — inA~^ , = d — inA~ 
and A'^ ,A~ are the (1,0) and (0, 1) parts of A. 
Next we define the following Laplacians: 

A = -(★!)★ D-h-D^-D*), I^^ = -2{-kD'^ i.D^ + i^D^i.), (94) 

which reduce to the Hodge de Rham and (anti)holomorphic Laplacians when acting on uncharged 
forms. The fact that these three Laplacians are equal when acting on uncharged forms is a well- 
known result for Kahler manifolds. Interestingly, this fact is also true when acting on charged 
two-forms. 

As discussed in the main text, for doubly transverse tensors, one may consider the hermitian 
and antihermitian cases separately. ^'^ In fact it turns out that hermitian doubly transverse tensor 
harmonics do not exist on CP^. This is easy to prove, at least in the uncharged case. Any hermitian 
tensor can be written as Yab = J^'^w^f)- Transversality implies Duj^ = 0, which reduces to du)^ = 
in the uncharged case. From the cohomology of CP^ it follows that uj^ must be proportional to J. 
But J is self-dual and u)^ anti-self dual hence = 0. 

For antihermitian modes = 0. Note that L^'^uj'^^ and Kg^u}~^ are the (0, 2) and (2, 0) parts of 
Yab^ and we can treat these separately. It is easy to show from the explicit form of L and K that 
transversality of Yah implies D^uj^ = 0. Let us focus on the (0, 2) part. Results for the (2, 0) part 
follow from complex conjugation. In jSHI the eigenvalue of such a two tensor was calculated for uo^ 
with charge n = 2, corresponding to Y^b of zero charge. We will generalise this to Yab of arbitrary 
charge n and thus uj~^ of charge n + 2 (since L has charge —2). We find that (dropping the subscript 
on henceforth) 

Ai{La'0Jcb) = Lac (A + 2 A) OJ,b (95) 

which can be proven by showing that 

A 2 4A , , 

/\U)ab = -'D idab + "^^afe (96) 

for an anti-self dual two form uj of arbitrary charge. The uncharged case follows from the Hodge 
de Rham Laplacian on two-forms. The charged case can be proven explicitly by showing that 
A — AhcIR = — + when acting on an anti self-dual charged two-form. For this computation 

^^Note that this is the only point at which we refer to the doubly transverse condition. Once we have divided 
into hermitian and anti-hermitian perturbations, the doubly tranverse condition is a consequence of the transverse 
condition. 
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one needs the explicit form of the Riemann tensor for a Kahler metric of constant holomorphic 
sectional curvature (which CP^ is an example of!), which can be found in ||26j. 

Following [2S| we define a differential operator E which when acting on a scalar gives an 
anti-self dual two form in the kernel of , just as we require. Explicitly, 

E(j) = D+i.D+{L(p). (97) 

This allows one to construct all such eigentensors of the Lichnerowicz operator from scalar eigen- 
functions of the Laplacian, which is the technique used in ^26,. In particular since we want uj to be 
of charge n + 2 we can always set uj = E(f)ri+A for some scalar field of charge n + 4. Next we need 
the eigenvalues of E(j)ri (since this is a two form we may choose any of the Laplacians). Using 

{D+,D-} = —J, *{JA*D+Lq))=-ii<D+Lc^ (98) 

one can show 

AE<f>n = [\n + 2(n - 1)A] Ecl>n (99) 

where —'D'^(j)n = Xni'n- Note this formula was proved in 2Q for n = 2,4. Putting these results 
together we find: 

Af (L>,6) = [A„+4 + 2A(n + 4)] L>efe (100) 

where uo = E(j)n+A- Thus to finish the computation of the eigenvalues of these tensor harmonics we 
need to know the eigenvalues of charged scalars. Fortunately this was done for CP^ in [2^, which 
uses the same conventions as we have. Thus first we discuss how to go between the conventions 
of j26j and ours. Firstly, we work with A = 6. Secondly our definition of charge is with respect to 
a potential A which satisfies J = ^dA. Thus the charge n in j26j is related to our charge m via 



n = The results of [21] give the eigenvalues of charged scalars: 



r + (|m| +2)A; + |m| (101) 



where A; = 0, 1, 2, . . .. We have introduced the modulus sign to take care of charges of negative sign, 
which can be deduced from those of positive sign by complex conjugation. Remembering to convert 
the charge we see that the eigenvalues of antihermitian (0, 2) tensor harmonics of charge m are given 
by: 

A = A„+6 + 8(m + 6) = /(/ + 4) - - 4m + 12 (102) 

where / = 2A; + |m + 6| and A; = 0,l,2,.... It turns out that / is the corresponding quantum number 
on S^. The eigenvalues for antihermitian (2, 0) tensors follow by complex conjugation and thus take 
the same form with m — > —m. 

Now we face a subtlety concerning the ranges of k. Obviously the eigentensors must be non- 
vanishing so we need to examine when E<j) = 0. Fortunately has addressed this issue and we 
can use the same technique. This involves introducing the adjoint of E, which we will denote by E^ , 
with respect to the Hodge inner product and thus maps two forms to scalars. An explicit expression 
for this adjoint is given in j26j : 

E'fuj = i.D-{K AirD-uj). (103) 
Note that Ecp = if and only if E'^ E(j) = 0. Using (|98|) one can prove: 

E^Ec^r^ = -(A„ + nA) [A, + 2(n - 1)A] (104) 

Note that this was proven in [22] only for n = 2, 4. Therefore Ecp^+A = if and only if: 

An+4 = -(n + 4)A, or A„+4 = -2(n + 3)A. (105) 
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Converting to our notation, we thus learn that one must exclude any values of k for which: 

= -4(m + 6), or A^+e = -4(2m + 9). (106) 



This is easily done and we find that one must exclude A; = 0, 1 for m < —6 and A; = for m = — 5. 
Pleasingly this implies / > 2 as one would expect for tensor harmonics. Once again, simply let 
m — > — m to deduce what happens in the holomorphic case. 

It is important to check that the modes we have considered are not pure gauge. We do this by 
computing the eigenvalues of pure gauge tensor modes with respect to A^. A pure gauge perturbation 
is of the form h^i, = V(^^j^) on the spacetime. Reducing to the base will give Yab = D^^^^f^y Using 
the gauge transverse and traceless conditions and working in the conventions of gives: 



(107) 



Now we specialise to antihermitian y^f,, so JY = ieY; this can be achieved by setting = ie(,. 
Then 

AiYab = A (^2 + e^^ Yab. (108) 

Converting to our conventions implies that pure gauge eigenvalues take the form: 

A = 2 (6 + em) . (109) 

By comparing these to (|102j) , and taking the ranges of k into account it is easy to show that the two 
sets of eigenvalues do not overlap so we do not need to worry about pure gauge modes. 

C Asymptotically flat case: WKB potential 

The maximum value of the function P defined by (|80|) is 
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where we have used uj < m^H and the upper bound on $7//. The expression in square brackets is 



positive and monotonically increasing (as a function of X) for X < 
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The expression on the RHS is monotonically decreasing for Y < y > from which we conclude 
Q > and hence Pmax ^ 
black hole is extremal). 



Q >0 and hence Pmax > 0. Equality occurs if, and only if, uj = mVlH and = \ / (jv^)^^ (i-^- the 
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